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Abstract 

We formulate a gradual dynamical transition from a power-law inflation era 
with a scalar field to a radiation era with no scalar field including inhomo- 
geneous perturbations to the Priedmann- Robertson- Walker universe. We 
show that for the cosmic microwave background radiation fluctuations this is 
excellently approximated by a sudden transition, with application of the Lich- 
nerowicz conditions, both for density and gravitational wave perturbations. 

I. INTRODUCTION 

Post-inflationary reheating in which the energy of a scalar field of the inflationary era is 
converted into the relativistic particles of the radiation era is of intrinsic interest, and there 
may prove to be tests of the mechanisms of inflation and reheating. A common scenario 
is that there is a period of inflation followed by reheating, consequently followed by the 
radiation era; density and graviton perturbations are supposed to be of quantum origin 
arising in the inflaton field(s) early in the inflationary era. It is an interesting question 
how the reheating affects the perturbations of various wavelengths, particularly the density 
perturbations. For those wavelengths A observable in the variations of the CMBR the phase 
At/A, where At is the reheat time, is extremely small. So the usual assumption is that for 
these wavelengths the effect of the transition on the perturbations can be well calculated by 
the sudden (instantaneous) approximation. In recent years there has been some discussion 
on this point JTJ-Q.It is the main purpose of this paper to support the usual assumption 
by illustration using a well known inflationary scenario followed by a particular reheating 
mechanism. 

The mechanism that we use is an old one [^], that of a friction term in the scalar field 
equation with the corresponding balancing term in the radiation fluid energy equation. We 
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do not postulate this as a close representation of the actual physical process of reheating 
or defrosting, which is anyway not yet determined, but use it as a stand-in providing a 
balanced transfer of energy-momentum from one phase to another in a mixed phase state. 
In this latter respect it is more physically motivated than reheating representations which 
just interpolate the scale factor of the FRW universe 

The new feature is that the calculation of the transition from a universe of all scalar 
field to one of all radiation fiuid is done explicitly for the perturbed as well as for the un- 
perturbed case. In the former case the synchronous gauge is used for calculation. We find 
that the numerical results for the resulting perturbations in the radiation era are the same, 
to many significant figures, as those got in the sudden approximation by application of the 
Lichnerowicz conditions 0] on a hypersurface of constant energy density, as formulated by 
Deruelle and Mukhanov 0. This result serves to confirm the hypothesis (and its imple- 
mentation) that the sudden phase transition happens on a hypersurface of constant energy 
density even where it is between entirely different types of physical states - in this case a 
scalar field and a standard type of perfect fluid. 

In the following section we outline the model, setting out the scalar potential. This 
is adjusted to give a period of power-law inflation, followed by a transition to zero scalar 
fleld, with the transition being mediated by a friction term. In section |T| we flrst derive 
the resulting equations for the reheating transition in which the scalar fleld matter of the 
universe is wholly converted into a radiation fluid and the universe enters the radiation era; 
secondly we derive the equations for the development of the density perturbations, arisen by 
earlier quantum fluctuations, through the gradual reheating transition. Section |^ is mainly 
on the numerical solution of these equations to give the resulting radiation era parameters; 
for the perturbations we treat only wavelengths important for the observed fluctations of the 
cosmic microwave background. We also derive the corresponding radiation era perturbation 
for the sudden transition approximation, so as to compare this with the result of the gradual 
transition. We show that the two results are the same, to extremely good accuracy. The 
computations are done in the synchronous gauge, but the comparison of results for the 
sudden and gradual transitions are done using gauge invariant variables . We also include 
a discussion of the effect of our gradual transition mechanism on non-physical so called 
'synchronous gauge modes' - eliminable by a coordinate transformation using the residual 
gauge freedom of synchronous gauges 0,51 • 

The metric tensor perturbations - gravitational waves - are a separate and much easier 
case and in [IV (J| we show analytically that the sudden transition approximation is valid also 
in this case. 



II. MODEL FOR INFLATION AND REHEATING 

We choose a potential for the scalar fleld (j) which gives constant power law inflation for 
large and which, in conjunction with a friction term, converts the scalar fleld energy into 
a perfect fluid as decreases. We take the fluid to be the usual one of the radiation era, so 
that the pressure and density are related by p = The potential is given by 

V = Uexp{X<f)),{X>0};<f)><f)A, (1) 
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V = yo[exp(-2i/0) - 2exp(-i/0) + 1]; (/>< 0^- (2) 

and is illustrated in the Figure for the particular case where Vq = l,v = \/2 and A = .45. 

The constants U, Vq, A, v are adjusted so that V, dV/dcf) are continuous at = 0^. The 
curvature, cPV/dcf)'^, changes sign at 0^ < 0a so that it is negative for (f) < (j)B where a 
friction term fsJdPvJd^, where / is an adjustable constant, is added to the equations as 
detailed below. 

The scalar field decreases with time and in the power law inflation region, > 0^, the 
cosmic scale factor 

a(T) oc (Tp - rf cx t« (3) 

where r is conformal and t cosmic time. Tp is a constant, A < -\/2 and 

p = 2/(A^-2),g = 2/A^. (4) 

For (j) < (pA the field accelerates until some (j) < (f)B where the friction becomes large enough 
to slow it down; finally the field energy is all converted into fluid energy. 

This model is just one realization of inflation followed by reheating and is designed to 
track coupled metric,scalar field and density perturbations through inflation to the end of 
reheating. These perturbations are supposed to arise from quantum fluctuations early in the 
inflationary era. They are accompanied by quantum fluctuations with graviton production 
giving rise to gravity waves which we shall likewise track through to the radiation era. 



III. INFLATION AND REHEATING EQUATIONS 
A. The unperturbed equations 

Wc write the equations for the case where the radiation fluid and the scalar fleld are 
coexisting; the pure scalar field era and the radiation era are just special cases of these 
equations. For 4> > 4>a the model has an analytic solution but otherwise a numerically 
computed solution is necessary, even for this unperturbed case. 

The Einstein equations are 

K = SttG (6) 

the right hand side being the sum of the energy-momentum tensors of the scalar field 
and of the fluid, with density p and pressure p = bp and 4- velocity u^^: 

t;. = M,. - g,A\9"^4^,o.<t>,p + vm (7) 

= P9ixv + (P + p)u^,u^ (8) 
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To solve these equations the scalar field equation is also requisite. To the usual scalar 
equation we add a friction term so that it becomes 

(fcpldr'^ + 2aHd(p/dT + aVi + aVfdcp/dr = (9) 

where Hubble, H = {da/dT)/a? and we have used the notations 

Vi = dV/d(j), (10) 



Vf = f^dW{(j))/d(j)^, (j)<<j)B,Vf = 0,<j)> (j)B. (11) 

From the Einstein equations we obtain 

H' = l[{^r/2a' + V + p]. (12) 

The Bianchi identities from the Einstein equations combined with the scalar equation yield 
the fiuid energy equation; 

dp/dr = -3aH{p + p) + a-^fid^/drf (13) 

Eqs. (|,|1^,0) are sufficient to solve the model; we note that as decreases kinetic energy 
is taken from the scalar field and given to the fiuid. The Einstein equations remain valid. 
For calculation it is convenient to use 

a; = Ina (14) 

as evolution parameter instead of r and we use the notation ' = d/dx. The equations of 
motion then become 

+ ^0' + 3(/)' + Vi/H'' + Vf(i)'/H = 0, (15) 

i/2 = «:(r + p)/(3- V), (16) 



p' = -Ap + HVf(j)'\ (17) 

where in the last equation we have used p = \p for the relativistic fiuid. In these equations 
V,Vi,Vf are given functions of cj) only, and the solution of the equations gives p{x), (f>{x), H[x) 
as functions of x = Ina only. We can find the value of r as a function of x and thus of a by 

r = |'^exp(-a;). (18) 

We note that the equations scale in the following sense: if we have a solution with a 
certain Vq, Eq.(^, with a certain 0(x), p(x), H{x) and c is any constant then c^Vq, 0, c^p, cH 
is also a solution. This is convenient because if we have a solution with a particular Vq then, 
by adjustment of c, we can get a solution appropriate to the era we wish to consider. 
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1. (j) > (pA'j power-law inflation 



In this region p = Vf = 0,V = U exp(A0), Vi = f/A exp(A0) and Eqs. ([T5|JT^) yield 

07(1 + 3(0' + A/«:)=0. (19) 

D 

The simplest solution is 

0' = --,i/2 = 2^7exp(A0)/(6-— ). (20) 

K K 

Letting A = A^/k it follows that H oc exp(— ^), r — Tp oc exp[(|- — l)x], a oc (r — Tp)P,p = 
2/(A — 2), Tp = constant. With A < 2 this is the well known power-law inflation due to an 
exponential potential. 



2. (l)<(j)A 

In this region we need to use numerical compution to solve even the unperturbed equa- 
tions whose solution we shall denote by 0o,Po- In our numerical computation the unper- 
turbed solution and the perturbation to it are evolved simultaneously, and we proceed now 
to the formalism for this latter. 



B. The perturbed equations 

Perturbations to the energy-momentum tensor density are coupled to scalar perturba- 
tions of the metric tensor. For these we work in a synchronous coordinate system and follow 
in the main the formalism of Grishchuk The metric is 

ds'^ = —a^dr'^ + a^[5ij + hij{x)]dx^dx^ . (21) 

hij{x) are the relevant perturbations which can be expressed as a linear superposition for 
the various wave numbers k of 

hij{k, x) = h{k, T)Q5ij + hi{k, T)k-^Q,^,j, (22) 

where Q is a superposition of the spatial wave function solutions exp(±zk.x) so that Q^ij = 
—kikjQ. The total perturbation can be found by solving for h, hi and this we shall do. We 
write the k-components of the perturbation to the energy-momentum tensor as 

T°o = Pig, r° = -T^o = a-2<Q,, (23) 

n, = {p,+pi)Q6i + k-%Ql. (24) 

where a = dx/dr = ^/a = aH, and where pi and pi,pi are density and pressure per- 
turbation components appropriate to k. Then the perturbed Einstein equations where the 
variables are functions of a; = In a are 
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3h' + [k/aYh — h[ = K{a/aY pi, 



(25) 



h' = Ki' (26) 

- h" - ^h' - 3h' = K{a/a)%, (27) 
H 

\[K + + 3/.; - [k/afh] = K{a/afp,. (28) 

In what follows we shall denote the k-component perturbation to the scalar field by 0i and 
the k-component fiuid perturbations by pi,Pi,^' and pi. For the perfect fiuid, which we 
assume, pi = and the corresponding quantity vanishes for the scalar field; so pi = 0. The 
right hand sides of Eqs.(p5|-p8D are given by: 



(a/«)2pi = + (0il^i + p,)/H^, (29) 
e = e- 0001, (30) 
ia/a)% = - (01 - ^pi)/H', (31) 

{a/afpi = 0. (32) 

Inserting these into the preceeding equations we have a set of differential equations capable 
of solution when we have added the following perturbed scalar field equation; 

0'i' + + 30; + + [^] + liSh' - hM = -[^0001 + Vf<Pi]/H (33) 

This equation combined with the (perturbed) Bianchi identities yields the perturbed fiuid 
energy equation as: 

p[ + 4pi - {k/an' + Ipoi^h' - h[) = [^0oVi + 2Vf<PMH (34) 

As in the unperturbed case the physical assumption lies in the addition of the friction 
term to the scalar equation; the Einstein equations are preserved leading to the perturbed 
fiuid energy equation, including the appropriate friction term which balances the one in the 
scalar equation. 
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IV. GRADUAL AND SUDDEN TRANSITIONS 



A. Inflationary era 

The initial conditions for solving the equations in the transition region, 0a > > 0, are 
given by the solution at the end, = 0a, of the power law inflation in the region > 0^. In 
that region we already have the unperturbed solution in |111 A 1| and shall now use well-known 
analytic methods to get the perturbed solution. 

We shall use dotted quantities to denote differentiation with respect to r; for example 

h = dh/dr (35) 

From Eqs.(g|,|0|) 

/€01 = -/i/0o (36) 
It can be shown that for the scalar field only case 

f, + fi[e-{a:/^)/{a^)]=0 (37) 

/i= ^-(/i + a7/i), (38) 
"\/7 

7 = l-d/a^ (39) 



where, as defined in |IIIB| , a = a/ a. The scale factor is given by Eq.(|^) so that 7 = 
(1 +p)/p = constant. 

Integration of Eq.(|38D gives h as 

h = -{Vl f\ik,T)dT + Q}, (40) 

where Cj is an integration constant subsuming, when 7 is constant, the lower limit of the 
integration, h = -Ci is a solution of = which can be eliminated by a coordinate 
transformation allowed by the residual gauge freedom of the synchronous gauge [|I| ; it gives 
zero contribution to gauge invariant variables We drop this term here but we shall 

revisit it later in conjunction with the results. 

All the above development was non-quantum mechanical. We shall now briefly remark 
on the quantum basis of the phenomena. We denote the corresponding quantum field theory 
quantities by a tilde: 

h = -^ r mr)dT (41) 
a J 

- — 3 [cfc/ii(fc, r) exp(2k.x) + h.c] (42) 
(27r)2 v2/c 
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where Ck is a quantum annihilation operator, [cfc, c^/] = 5^(k — k'), and fii{y),y = |/c(rp — r)|, 
is that solution of the Bessel equation ( |57| ) such that 

^iM^e-'y,y^oo (43) 

thus corresponding for large k to the usual mode function of quantum mechanical plane 
waves. 

is a normalization factor whose determination gives the absolute magnitude of the 
observed density perturbations (in terms of the model parameters) given the assumption 
that these come from a primordial vacuum with zero quantum occupation number. N has 
been determined as 



N = V2i^ = VI671G. (44) 

For this paper, which is mainly about sudden transition versus gradual transition am- 
plitudes, we do not need to discuss the quantum aspects further. 

We can now proceed with the evaluation of the metric scalar components at the end of 
inflation. The exact expression for the solution /ii is, with n = ^ — p, 

pTvy 1 1 

^J'l{y) = y Y^'^"" ~ ^^""^ exp[-i(-n7r + -tt)] (45) 

Exact power-law inflation ends and the transition region begins at = </)^; other variables 
at that surface, A , we denote by the suffix A. For example: yA = \k{Tp — ta)\- For small 
enough values of yA we can expand /z given by Eq. (^5D in an ascending power series and 
just consider the first one or two terms. Our first interest is in values of k of importance in 
COBE and other observations. These are of the order of magnitude (the relevant meaning 
here is being within a few factors of ten) of 

ki = a\/ai = aiHi (46) 

where ai,Hi denote the values of the scale factor and Hubble at the time, ri, when the 
matter era begins. Thus for such values of k the values of yA are of order 10^", n > 10. 
Then we expand /x in a power series with leading terms 

^'^(y) = M(p)iy"~ 2i2pli) ^^^^ 

M{p) = -2-Pr(i - p) exp(-z(l - p)/2)/V^. (48) 
The corresponding expansion of h from Eq. (^0]) , with Ci = 0, is 

h = (^a)-^M(p)[yP - ^ ] (49) 

VV I J WLi/ 2(2p+ l)(p + 3) ^ ^ ^ 

We need these expressions at y = yA to find the initial values for the numerical solution of 
the coupled differential equations in h, /i^, 0i, pi. Noting that the first term in the expansion 
of /i 
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S = {^a)-^M{p)y^ (50) 
is a constant, we find to first order in y\ the initial values 

h = T. + ylPj: (51) 

h' = ^ylP^E (52) 
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n = ^ylPi^ (54) 



^; = -01 



where 



(55) 

p'l = (56) 
Pi = -(p+l)/[2(p + 3)(2p+l)] (57) 
B. Reheating and the radiation era 



In Eqs. (^5| - p8D h occurs only when multiplied by {k/a)"^ = y^/p^. This and the initial 
conditions above lead us to rewrite the equations in terms of renormalized variables which 
are of order zero in y'^. These are: 

h = {aA/k)\h - S)/S , h[ = {a/kfh'J^ , 0i = {aA/kf<P/^ , pi = {aA/kfp/^ , 
6 = {aA/kniT.. 

Except for h'l the renormalization is by the same constant factor coinciding with the 
factor for h'l at A. The initial conditions for the hatted variables are given by the replacement 
y\Ti — > p^ on the right hand sides of Eqs. (|5T| - |56|) . The hatted versions of the last two Einstein 
equations are 

h" + + 3)h' = -k[<P',4>[ - (01^1 - Ipi)/H% (58) 
CKy + {l-^)K-l = (59) 
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We shall supplement Eqs . (1581 , 1591) by Eqs. ( pOlplD below, which are the hatted versions of 
Eqs.(33|.34). We shall omit the explicit k"^ terms because these are of order {k/aY = y'^/p'^ 
smaller than the other terms in the equation for both cases. This is a consistent, controlled 
approximation, extremely good for COBE and related observations. 



+ (| + 3)0; + |l 



^ + l[3h'-{{aHU/{aH)rh[] 



4Vf 



+ Vj^',yH (60) 



p\ + 4pi + ^po[3/^' - {{aH)j,/{aH)fh[] = [^cj^'^k + '2Vfct>'ok]H (61) 



Eqs.(58)-(|6l|), with the initial conditions specified above, have to be solved for the pertur- 
bations. It is necessary to know the unperturbed solutions 0o, Po and H which enter the 
equations; these are given by eqs. (|l5|) - ([17D and their initial conditions at the beginning of 
reheating are given in [111 A 1| . We evolve all the variables, both the perturbed and the unper- 



turbed, together as a set of coupled differential equations with evolution parameter a; = Ina. 
The end of the transition period and the beginning of the radiation era are signalled by 0o 
arriving at a negligable value; 0i accompanies 0o into oblivion. 

Our objective is to compare the radiation era amplitudes resulting from this gradual 
transition to those resulting from a sudden transition. 



1. Sudden transition to the radiation era 

For this we use the matching conditions given by Deruelle and Mukhanov ; conditions 
(their eqs. 4.5a,b) in the longitudinal gauge can equivalently be given in terms of the gauge 
invariant variable $. They are that $ and F are continuous across the matching surface 
r = Ta: 

$+ = <l>(r2+) = <l>(r2_) = $^ (62) 



F+ = F(r2+) = F(r2-) = F_ (63) 

where 

F = ($' + $ + e^$)/7, (64) 

e = k/{aV^). (65) 

Additionally, more obviously, a and a = a/ a = aH must also be continuous across the 
surface. 

Density perturbations can be described by the mode function i/, corresponding to the /ii 
of the inflation era, such that 

j) + ifcV = 0. (66) 
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In terms of the synchronous gauge variables 

u = -{h + a-fh) = a{h' + -fh) (67) 
a 

and the gauge invariant variable $ is given by ^ 

$ = - {a/kYh'i (68) 

so that from Eqs.( ^^ - |28| ) 

(^ = {u-u')/ae^ (69) 

Non-physical synchronous gauge modes (that is of the form h = a/a x constant) give zero 
contribution to both u and $. We write the solution of Eg. as 



so that 



Ccos(A;(r - T2)/V3) - ^sin(A;(r - r2)/V3). (70) 



^+ = {e2S + C)/ia24), (71) 

r+ = -[Cil - el) + e2Sil - lel)]/ia2el). (72) 

From the continuity of $ and F across the interface these expressions enable the determina- 
tion of C, S and thus the density perturbation in the radiation era. So we need the values 
of $_ and F_ from the inflationary era solutions. From Egs.(|51|, p3| , |68D , interpreted now as 
being at the surface r = T2, we find that to lowest order in (k/a)^ 

$_ = ^^S; F_ = -^—J:. (73) 
2p + 1 ' 2p+l ^ ^ 

Solving Egs.(^,^) straightforwardly gives to lowest order 

C = -eaS = 2a2S. (74) 

For k relevant to CMBR fluctuation observations €2 is of order 10~" where n > 10 and thus 
away from the interface the term in S dominates the radiation era density perturbations. 



2. Gradual reheating transition 

For our reheating model, or any such continuous transition, no guestion of sophisticated 
implementation of Lichnerowicz conditions at an interface arises. There is no interface; the 
physical variables arising from the solutions of the eguations of motion are all continuous 
through space-time. 

We shall consider the case where the transition ends at the same surface, r = r2, as that 
of [IV B 11 (f) is negligable for t > T2 and the system is in the radiation era. The transition 
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begins at the end of pure power law inflation, r = ta- We note here that S of Eq.(0) is 
a constant of the power-law inflation era. Below we shall refer to a S evaluated at at 
the end of power-law inflation at the beginning of the gradual transition; while in [IV B 1 



we 



referred to S evaluated at that appropriate (sudden) end of power-law inflation, T2; these 
have the same name and are indeed the same numbers. 

The radiation era analysis is of course precisely the same as in |IV B 1| , with Eqs.( |68| , |69D 
giving 



3$ + $' 
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so at the end of the transition, t = T2, which is the beginning of the radiation era 

C/a2 = 3<l>2 + $2, 



(75) 



(76) 



(C + e25)/(a2e^) = $2. (77) 

Putting into the RHS of the last two equations the computed values at the end of the 
transition we can find C and S. In fact, because of the very small magnitude of €2, the last 
equation implies that to lowest order in €2 

628 = ~C (78) 

since otherwise $2 would take a large value unattainable in the transition. The result we 
have arrived at is that for the continuous case the density perturbation in the radiation era 
is determined by the value of 3<l'2 + '^'2 cit the end of the gradual transition. 

By numerical integration of Eqs. (|5^ - |61|) with the stated initial conditions we determine 
$2 and $2 terms of S using Eq . (|68|) . We can then compare the results for 5* and C, 
determining the radiation era perturbations, from the sudden and gradual transitions. 

3. Comparison of sudden with gradual transition 

For the sudden transition we have, Eq. (|7^ , 

- e25/a2S = C/a2S = 2 (79) 
For the gradual transition we denote 

result = (3<l>2 + $2)/^ (80) 

so that 

— e2S/a2T, = C/a2T, = result (81) 
Thus if result = 2 the two transitions agree. 

We have calculated the gradual transition for two values of the inflation power-law ex- 
ponent, p, each for a number of values of the friction coefficient, /, with the results shown 
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in the Tables. It may be noted that for the larger values of / ,(/ > 1) quoted the system 
grinds to a halt at the bottom of the potential well with no significant oscillations; this 



of course does not resemble any usually envisaged reheat mechanisms For smaller / 
oscillations occur and it is found, as we would expect, that their number increases rapidly 
as / decreases. Our numerical computation methods fail for / significantly smaller than the 
last value in the Tables, / = .01. 

As stated previously the transition is deemed to have come to an end when the value of 
00 is sufficiently small. There is not a precise criterion and we display results for various 
values of (pQ^end). We see that always result = 2 to very high accuracy, in agreement with 
the sudden transition. 

This establishes the agreement (for appropriate wavelengths) in the following sense: 
whether there is a sudden transition at red-shift Z2 or our type of gradual transition ending 
at Z2, the resulting radiation era scalar perturbation is the same. 

4- Synchronous gauge mode 

As set out in Eqs.(^)-(|40|) there is in the inflationary era (as in other eras) a 'synchronous 
gauge mode' given by 

h = -C, (82) 

where Cj is a constant. This mode is non-physical as it contributes zero to the gauge 
invariant variable $, given by Eq.(^), as well as to other gauge invariant variables, and 
also it can be eliminated by a coordinate transformation making use of the residual gauge 
freedom of synchronous gauges. 

Whole these remarks justify our ommission of this term, we have investigated further so 
as to provide a comment on the mechanism of energy transfer inthe gradual transition of 
our model. Our theory is a linear perturbation theory , so we can completely investigate the 
consequences of a synchronous gauge mode separately from the physical modes previously 
considered by inserting solely Eg . (^2)) and its derivative and 

= (^kjaf-Ci (83) 

(X 

as the initial conditions in Eqs.(|58|)-(|6lD. These correspond to zero $ in the inflation era. For 
a perfect, that is coordinate independent, perturbed transition mechanism in Eqs.(|5U[) and 
(|6T|) the resulting $ in the radiation era would be zero. To give an idea of the magnitudes 
involved we can put, formally, (^)yiCj = S at the beginning of the transition - that is as if 
this quantity began the transition with the same magnitude as the physical mode previously 
considered. We illustrate by giving two results, corresponding to the first and last lines of 
Table I, for the value of $ at the end of the transition: 

/ = 4^<l>2 = Sx 10-^ (84) 
/ = .01 ^ $2 = 4.36S X 10"^ (85) 
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The results for the physical mode were 2S for whatever /. Eqs. (|8^ ) and (|^) illustrate 
firstly that there is a coordinate dependence in the transition mechanism when applied to 
the non-physical, coordinate dependent, perturbative mode and secondly that it decreases 
with /. 

It is probable that other transition mechanisms expressible as terms in the Lagrangian, 
such as some of those which have been invoked to investigate particle production in the 
reheat phase [^, would transmit only coordinate independent information into coordinate 
independent variables. However investigation of these would involve different formalisms and 
calculations. We would not expect any different answers for the physical modes. 



C. Gravitational waves 



Tensor perturbations are very much simpler to treat than scalar perturbations. After 
their generation by quantum fluctuations their development is governed solely by the cosmic 
scale factor a. Corresponding to Eq.(^) we have in the power-law inflationary era 



2 ^3^ 

hij = V1QttGJ2 / ,^ ,3 /— , M xkefj{k)^i{k,T)exp{'ik.:)^) + h.c], (86) 



t^J (27r)f v^a(r) 

where a\k is the annihilation operator for the graviton with polarization A and wave number 
k, and the polarization tensor satisfies £ij{k)efj{k) = 25\\'. /ii given by Eq . (|47|) is a 
solution, appropriate to the quantum mechanics as in Eq. (^3|) , of the tensor mode equation 

'fi + /i(A;^ - a/ a) = 0. (87) 

This holds for any cosmic era, whatever may be the dynamics responsible for the particular 
form of a(r) and in the radiation era, beginning at a = 02, which defines r = r2, we can 
write its solution as 

fi = G+ cos{k{T - T2)) - G- sin(A;(r - rs)), (88) 
where (7+ and are constants to be determined by continuity, as set out below. 



1. Sudden transition 

If we assume a sudden transition from inflation to the radiation era at a = 02 the 
matching conditions that determine G+ and G- are that fi and /i (equivalently fi') must be 
continuous across the interface; the continuity of a and a is also a condition but we already 
hold these enforced as in PTB| 11,0. Up to the second term of the power series expansion 



H is given by Eq. (|47|) and /x' by 



and the required continuity gives 
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= M(,mi - (90) 

Due to the extreme smallness of ?/2 for waves of observational interest dominates and to 
lowest order 

G+ = 0, (92) 

G.=pM{p)yP,-\ (93) 



2. Comparison of gradual with sudden transition 
At the end of power-law inflation, beginning the gradual transition 

I^M = M(p)All - (94) 

,'(r.) = M(p).5|l-|±|||. (95) 

Also we shall below use the fact (see Eq. (pOD ) that 

Sg == M{p)yP/a = ^T. (96) 

is a constant of the inflationary era. It is convenient to consider the function X = fi/ a so 
that from Eq.(|^ 

A + 2Xd/a + P\ = (97) 

Generally the solution is not trivial, but working to lowest order, putting k"^ = 0, we get 

a^X = constant (98) 

and, again to lowest order, at the end of the power-law inflation 

A = M{p)yP/a ^ Aa = 0. (99) 

Then A is constant throughout the transition, /i oc a, as it also is in the inflationary phase, 
though for different a(r); in both cases A = Eg-. This implies that, to lowest order, where fi 
and fii are the gradual transition and inflationary mode functions respectively: 

/i(a2) = /i'(a2) = /ii(a2) = /i'i(a2) (100) 

since a' = a. Thus Eqs.( p2| , |93|) hold as in the sudden transition case. 

This establishes the same result for tensor perturbations as was established for scalar 
perturbations, and in the sense stated at the end of section [IV B 3| . For scalar perturbations 



it was likewise justiflable to neglect terms of order {k/aY in certain equations for waves with 
k relevant to current observations. We may note that this easy way of flnding the result in 
the tensor case is not only due to the relative simplicity of the governing equation (0) but 
also to the assumption that the inflation is power-law. However we would expect it to hold 
more generally. 
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V. DISCUSSION 



As the transition proceeds both the unperturbed, 0o, and the perturbed, 0i, scalar fields 
diminish so that at some sufficiently small value we can use the perfect ffuid model of the 
radiation era. In the latter part of the radiation era, as the the matter era is approached, the 
amplitude of the density perturbation is specified by the single quantity result = 3$ + $' = 
3$ + {aH)~^d^/dT, evaluated at the end of the transition. 

An extension of the length of the transition, specified in the Tables by the number 
of e-foldings, can always give values of |0o| as small as we please until the computation 
runs out of numerical accuracy. We have chosen the varying examples in the tables to 
demonstrate that the lack of precision in the |0o| criterion is not important in the sense that 
result = {3^2 + ^2)/^ very near to 2 for many different 'end' values of |0o|- 

Though $2/^ and $2/^ vary in the third decimal place the result is much more constant, 
and we can understand this: In the radiation era by Eqs. ( [75| , [7DD (3$ + <I>')/S = —eS/aT, + 
O(e^) and the first term is constant since a (x r and e oc r and an equivalent result would 
apply towards the end of the transition region. Thus for the above calculations with small k 
the value of 3$ + $' is very constant over many e-foldings of the cosmic scale factor a after 
the scalar field has decreased in magnitude by more than about 4 powers of 10. This is a 
property of the early part of the radiation era and not a property of the particular transition 
mechanism adopted. 

Having concurred with some other authors |0-^,|T3[ , through different considerations, that 
transition period detail has no inffuence on density perturbations important for the CMBR 
ffuctations, we may ask for what wave numbers k the detail is going to be important. The 
answer is likely to vary considerably with the assumed physics of the early universe. So we 
can only give an answer for power-law inflation followed by reheating. For illustration it is 
simplest to consider the case of gravitational waves as in |1V C[ To get the gradual-sudden 
equivalence there required that ]?k^ jo? <^ 1 throughout the transition region. We can 
conclude that for k/ki ^ 10^^ the sudden result is still reasonably good; here h is the 
Hubble parameter, fci, given by Eq.(46, for redshift zi = 10^'^, being the value around which 
we did our tabulated calculations above, specifies the region of wave numbers of importance 
in CMBR fluctuations. Translated into energy terms this means that gravitational waves 
with k < 10~^^h'^eV would not be affected by reheating details. A calculation with density 
perturbations gives similar results. 
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TABLES 



TABLE I. result = (3$2 + ^2)/^ the 'end' (see text) of the gradual transition for values of 

the friction coefficient, /; R^f) is the ratio of the unperturbed scalar field value, (f)o, at the transition 
'end' to that at the beginning; x{end) gives the number of e-foldings of the transition. Conformal 
time inflation power p = —1.113. 



/ 


x{end) 


\R<p\ 


$2/S 




{result — 2| 


4.0 


3.2 


0.29 X 10"^° 


0.66162 


0.01515 


<10-8 


2.0 


3.0 


0.20 X 10-23 


0.66262 


0.01214 


<10-8 


1.0 


3.0 


0..30 X 10-^^ 


0.66343 


0.00971 


<io-8 


0.5 


3.0 


0.34 X 10"^° 


0.66337 


0.00990 


< 10-8 


0.1 


3.5 


0.18 X 10-^ 


0.66450 


0.00650 


1.1 X 10-8 


0.05 


4.0 


0.14 X 10-6 


0.66534 


0.00399 


<10-8 


0.01 


5.0 


0.54 X 10-^ 


0.66540 


0.00380 


<10-8 



TABLE II. result = (3<I>2 + *I'2)/^ t^*^ 'end' (see text) of the gradual transition for values of 
the friction coefficient, /; ii^ is the ratio of the unperturbed scalar field value, (po, at the transition 
'end' to that at the beginning; x{end) gives the number of e-foldings of the transition. Conformal 

time- iufialion power p = —1.05. 



/ 


x{end) 




$2/S 


$'2/5] 


\result — 2| 


4.0 


4.0 


0.61 X 10-^° 


0.66236 


0.01293 


<io-8 


2.0 


4.0 


0.21 X 10-1^ 


0.66476 


0.00572 


< 10-8 


1.0 


3.0 


0.26 X 10-^ 


0.63595 


0.09215 


2.5 X 10-8 


0.5 


4.0 


0.11 X 10-^^ 


0.66511 


0.00467 


< 10-8 


0.1 


4.5 


0.50 X 10-9 


0.66566 


0.00303 


<10-8 


0.05 


5.0 


0.48 X 10-^° 


0.66605 


0.00185 


<10-8 


0.01 


6.0 


0.15 X 10-^° 


0.66608 


0.00176 


<10-8 
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